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Abstract 

> 

\f} | After some generalities on neutrino oscillations and on neutrinos, 

■ the recent experimental results presented by Soudan 2, MACRO and 

SuperKamiokande at the Neutrino'98 conference are summarized and 



O | discussed. 

ON 
On 



1. Introduction 



45 

Atmospheric neutrinos are produced in cosmic ray interactions in the 
^ ■ upper atmosphere: a high energy primary cosmic ray, either proton or nu- 

cleus, interacts producing a large number of hadrons, in particular pions 
and kaons. These can decay giving rise to muons and muon neutrinos; 
also the muons decay yielding muon and electron neutrinos. In this sim- 
plified picture, the ratio of the numbers of muon to electron neutrinos is 2, 
(N^+N Vti )/(N Ue + N Ve ) ~ 2, and (N v /N v ) ~ 1, see Fig. 0. One may consider 
that these neutrinos are produced in a spherical surface at about few tens of 
km above ground, and that they proceed at high velocity towards the center 
of the earth. 

The neutrino flux has been computed by a number of authors. At low 
energies E v ~ 1 GeV, the absolute numbers of predicted neutrinos differ 
by about 30%@; at higher energies, E v > 10 GeV, the predictions are more 
reliable, with a systematic uncertainty of about 15%, that is about 1/2 of 
that at low energies!. However the predicted relative rates of and v e and 
the shape of the zenith distribution have a considerably lower systematic 
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Figure 1: (a) Illustration of the production, travel and interactions of atmospheric 
muon neutrinos; (b) interaction of a primary cosmic ray, production of pions and 
their decays leading to the atmospheric v e i v u- 



error. Other sources of systematic uncertainties are due to the knowledge 
of the neutrino-nucleon cross sections and to the propagation of muons and 
electrons in different materials. 

Since the 1980's several large underground detectors, mainly designed to 
search for proton decays, have studied atmospheric neutrinos. These detec- 
tors are located below 1-2 km of rock and they may detect neutrinos coming 
from all directions. Via charged current (CC) interaction the gives rise to 

and thus to a track, the v e to e~ and thus to an electromagnetic shower. 
The produced hadronic system is not observed in most cases. 

The detectors can be roughly classified as tracking calorimeters and as 
water Cherenkov detectors. The water Cherenkov detectors Kamiokandei 
and 1Mb! reported anomalies in the ratio of muon to electron neutrinos, while 
the tracking calorimeters NUSEXI and FrejuaO did not find any deviation. 
Also the Baksani scintillator detector did not see any deviation. 

Later Soudan 2! and MACRO0 reported deviations; in 1997 the neutrino 
anomaly was also reported by SuperKamiokande0. 

At the Session of June 5, 1998 at the Neutrino'98 Conference in Takayama, 
Japan, new, higher statistics data have been presented by the Soudan 20, 
MACROS and SuperKamiokandeB collaborations. The experiments 



ts con- 



2 



firmed the neutrino anomaly and opened a wide discussion on the subject. 

The main purpose of this lecture is to review and discuss the presentations 
at the Neutrino'98 Conference. The Sup erKamiokande results have been 
presented at this School by M. KoshibaEl Some results have been updated 
at the 1998 HEP Conference in Vancouver, CanadaB. 



2. Neutrino oscillations 



If neutrinos have non-zero masses one has to consider v e , v T and their 
combinations. The u e , v T are the appropriate particles to consider in weak 
decays, for example 7r + — > / u + + z/ At , and in charged current (CC) interactions 
which lead to their detection, e.g. + n — > ^ + p; in technical terms one 
says that they are weak flavour eigenstates. Instead in the propagation in 
vacuum the appropriate particles are the mass eigenstates Ui, z^, z/ 3 . The 
weak flavour eigenstates vi are linear combinations of the mass eigenstates 

3 

"1= U lm V ra (1) 

m=l 

In the simplest case of only two neutrinos (z/^, v T ) which oscillate with two 
mass eigenstates (z/ 2 , z/ 3 ) one may write 

= v 2 cos 9^ T + z/ 3 sin 6^ T ,^ 
v T = — z/ 2 sin 9^ T + z/ 3 cos 6^ T 

where Q^ T is the mixing angle. 

If the mixing angles are small one would have v e ~ u\, z/ M ~ z/ 2 , v T ~ z/ 3 ; 
in this case one may speak of the mass of z/ e which is about equal to that 
of Ui, etc. In the limit of zero masses the neutrinos become equal, and one 
does not need to introduce the z/ 1; z/ 2 , z/ 3 . If the mixing angles are large, the 
situation is different and one has to consider well separated the eigenstates of 
flavour from those of mass. In particular it would be inappropriate to speak 
of mass for the weak eigenstates v ei z/^, v T . 

In the case of only two types of neutrinos, and z/ T , one may easily 
compute the following expression for the survival probability of a beam: 

, . In It ^2 — E\ . In 2/1 . 27 A?7l 2 ■ L . 

P{Vn -> Vjx) = I- sin 29^ sin ( - 1) = 1-sm 26> Mr sin ( — ) 

(3) 



3 



The probability for the initial to have oscillated into a v T is: 

9 „ n . 1.27Am 2 • L. 
^r) = 1 - -> ^) = sin 2 2^ T sm 2 { ) (4) 

The mixing angle 9^ T and the mass difference Am 2 = Am^ may be deter- 
mined from the variation of P{yp — > z^ M ) as a function of the zenith angle G, 
or the path L, or the energy E u or L/E v . 



3. Early experiments 

The ring imaging water Cherenkov detectors Kamiokande and 1MB mea- 
sured Ufj, and v e CC interactions and found that the ratio of muons to 
electrons was smaller than expected! In the Kamiokande experiment, neu- 
trino interactions were classified using the shape of the Cherenkov rings on 
the phototubes on the wall of the cylindrical water container, and through 
the recognition of muon decays. The results were expressed in terms of 
the ratio R = Robs/Rue between R obs = (— ) & s of measured CC interac- 

' v Ve 

tion events and Rmc — (-)mc from Monte Carlo simulations. The sin- 
gle ratio {v^obsl (y^) MC , may be affected by large theoretical and system- 
atic uncertainties; in the double ratio most systematic uncertainties cancel. 
The measured double ratios from Kamiokande! and IMbU are RKamioka = 
0.60 ± 0.06 st at+sys and Rimb = 0.60 ± 0.06 stat ± 0.07 sys , respectively. The 
NUSEXQ and FrejusQ tracking calorimeter detectors reported for contained 
and semicontained events R Q b s ~ Rmc within errors. The measured double 
ratios are: R Nuse x = 1.0 ± 0.3, R Fre jus = 0.99 ± 0.13 stat ± 0.08 sys . 

The Baksan scintillation telescope detected a sizable number of upthrough- 
going muons arising from interactions in the rock below the detector!. The 
average energy for these events is considerably larger (50 — 100 GeV) than 
for the contained events measured by the other detectors. They did not find 
deviations from the predictions in the total number of events, but they find 
an anomalous angular distribution. 

Later, the Soudan 2 tracking and shower calorimeter confirmed an anomaly 
in the v^j v e ratio for contained events!. 

The MACRO0 collaboration reported in 1995 a first measurement of 
upthroughgoing muons coming from of E u ~ 100 GeV in which there 
was a deficit in the total number of observed upgoing muons, in particular 



4 





Gold Data 


Rock data 


Background 


Events 


MC 


Track 


95 


278 


18.5 ±6.1 


76.9 ± 10.8 


137.4 


Showers 


151 


473 


33.3 ±12.8 


116.3 ± 12.8 


133.8 



Table 1: Event summary for the Soudan 2 data. Track events are due to while 
showers to u e . The MC predictions were obtained using the v Bartol flux. 



around the vertical, and also reported an anomalous zenith angle distribution. 
The deficit was confirmed at the 1996 and 1997 conferences^. 

At the 1997 conferences, the SuperKamiokande Collaboration confirmed 
the Kamiokande and 1MB anomalies in the /i/e ratio for contained events^. 



4. Contained events and Soudan 2 

The Soudan 2 data were presented at Neutrino'98 by E. Petersonlll; im- 
proved data were presented at HEP'98 by H. Gallagher§. The Soudan 2 
experiment uses a modular fine grained tracking and showering calorimeter 
of 963 t. It is located 2100 m.w.e. underground in the Soudan Gold mine in 
Minnesota, USA. Its overall dimensions are 8m x 16m x 5m. The Central 
Detector is made of 224 calorimetric modules of lm x lm x 2.5m. The bulk 
of the mass consists of 1.6-mm-thick corrugated steel sheets which are inter- 
leaved with drift tubes. The detector is surrounded by an anticoincidence. 
This Active Shield detector covers the walls of the Soudan 2 cavern enclosing 
the Central Detector as hermetically as possible. 

The neutrino contained events were selected by a combination of a two- 
stage software filter and a two-stage physicist scan. The software filter rejects 
non-contained events by requiring that (i) no part of the event is within 20 
cm of the detector surface and (ii) no track is located or oriented in such a 
way that it could enter undetected in the calorimeter from a crack between 
modules. The last stage is the scan of real and simulated events. The MC 
simulated events were mixed with the real ones. All events were classified 
into one of three topologies: tracks (as due to CC interactions); showers 
[y e CC interactions) and multiprongs (interactions of all neutrino flavors and 
NC). The multiprong events are not considered at present. Finally, events 
without any hits in the Active Shield detector are called Gold Events, while 
events with two or more hits in the shield are called Rock Events. 
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Figure 2: Soudan 2 (preliminary) results on the number of observed events as 
function of L/E: (a) for Vn and (b) v e . Only statistical errors are shown; the 
dashed lines are the MC predictions using the Bartol flux. 



The data from a 3.89 kt-yr exposure are summarized in Table [I], together 
with the Monte Carlo predictions using the Bartol neutrino flux0. The 
Soudan 2 double ratio is R = (NjN e ) data / (NjN e ) MC = 0.64 ± 0.11±°;°5 
which is consistent with muon neutrino oscillations. 

For a smaller, high resolution sample, they are able to estimate the L/E 
for each event (Fig. |[), together with the MC predictions. Within limited 
statistics (after corrections, 60.8 CC events and 106.4 v e CC events), the 
preliminary data are consistent with an anomaly in the muon data and not 
in the electron data. For this data set R = 0.52 ± 0.09 stat . They add that 
Am 2 < 10 -3 eV 2 appears unlikely. 



5. Upward-going muons and MACRO 

The MACRO data were presented at Neutrino'98 by F.RongaEl. The 
MACRO detector is located in Hall B of the Gran Sasso Laboratory, with 
a minimum rock overburden of 3150 hg/cm 2 . It is a large rectangular box, 
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Figure 3: Sketch of different event topologies induced by interactions in or 
around MACRO. The stars represent scintillator hits. The time-of- flight of the 
particle can be measured for Up Semicontained and Up throughgoing events. 

76.6 m x 12 m x 9.3 m, divided longitudinally in six similar supermodules 
and vertically in a lower part (4.8 m high) and an upper part (4.5 m high). 
The detection elements are planes of streamer tubes for tracking and liquid 
scintillation counters for fast timing. The lower half of the detector is filled 
with trays of crushed rock absorbers alternating with streamer tube planes; 
the upper part is open and contains the electronics. There are 10 horizontal 
planes in the bottom half of the detector, and 4 planes on the top with wire 
and 27° stereo strip readouts. Six vertical planes of streamer tubes cover 
each side of the detector. The scintillator system consists of three layers of 
horizontal counters, and of one vertical layer along the sides of the detector. 
The time (position) resolution for muons in a scintillation counter is about 
500 ps (11 cm). Fig. |3| shows a sketch of the three different topologies of 
neutrino events analyzed until now: up throughgoing muons, semicontained 
upgoing muons and up stopping muons+semicontained downgoing muons. 

The up throughgoing muons come from interactions in the rock below 
the detector, with E v ~ 100 GeV. The muons (E^ > 1 GeV) cross the whole 
detector. The time information provided by scintillation counters allows the 
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Events Predictions (Bartol neutrino flux) 



detected No Oscillations With oscillations 

Up Through 451 612 ± 10A th ± 37 sys 431 ± 73 th ± 26 syi 

Internal Up 85 144 ± 36 th ± 14 sys 83 ± 21«, ± 8 SJ/S 

In Down+Up Stop 120 159 ± A0 th ± I6 sys 123 ± 31 th ± 12 sy , 



Table 2: Event summary for the MACRO analysis. The predictions with oscilla- 
tions are for maximum mixing and Am? = 0.0025 eV 2 . 



determination of the direction by the time-of-flight (T.o.F.) method. The 
data presented correspond to ~ 3.5 live years. 

The semicontained upgoing muons come from interactions inside the 
lower apparatus. Since two scintillation counters are intercepted, the T.o.F. 
is applied to identify the upward going muons. The average parent neutrino 
energy for these events is ~ 4 GeV . If the atmospheric neutrino anomalies are 
the results of oscillations with maximum mixing and Am 2 between 10 -3 
and 10 -2 eV 2 one expects a reduction of about a factor of two in the flux of 
these events, without any distortion in the shape of the angular distribution. 

The up stopping muons are due to external interactions yielding up- 
going muon tracks stopping in the detector; the semicontained downgoing 
muons are due to induced downgoing tracks with vertex in the lower 
MACRO. The events are found by means of topological criteria; the lack 
of time information prevents to distinguish the two sub samples. An al- 
most equal number of up stopping and semicontained downgoing events is 
expected, and the average neutrino energy for these events is around 4 GeV . 
In case of oscillations with the quoted parameters, a similar reduction in the 
flux of the up stopping events as the semicontained upgoing muons is ex- 
pected. No reduction is instead expected for the semicontained downgoing 
events (from neutrinos having path lengths of ~ 20 km). 

The data were compared with Monte Carlo simulations. In the upgoing 
muon simulation the neutrino flux computed by the Bartol group0 is used. 
The cross sections for the neutrino interactions have been calculated using 
the Morfin and TungH parton distribution set SI. The propagation of muons 
to the detector has been done using the energy loss calculation by Lohmann 
et aZ.Ei in standard rock. The total systematic uncertainty from neutrino 
flux, cross section and muon propagation on the expected flux of muons is 
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Figure 4: MACRO data, (a) Flux of the up throughgoing muons with > 1 
GeV vs. zenith angle O. The solid line is the expectation for no oscillations and 
the shaded region is a 17% scale uncertainty. The dashed line is the prediction for 
an oscillated flux with maximum mixing and Am 2 = 0.0025 eV 2 . (b) Probability 
contours for i/„ — > v T oscillations. The dashed lines are 90% and 99% CL contours 
calculated according to 21 . The best probability is 17%; iso-probability contours 
are shown for 10% and 1% of this value (i.e. 1.7% and 0.17%). 



~17%. Fig. |]a shows the zenith angle distribution of the measured flux of 
up throughgoing muons with energy greater than 1 GeV; the Monte Carlo 
expectation for no oscillations is shown as a solid line, and for a — > v T 
oscillated flux with sin 2 Id = 1 and Am 2 = 0.0025 eV 2 is shown by the 
dashed line. The systematic uncertainty on the up throughgoing muons 
flux is mainly a scale error that doesn't change the shape of the angular 
distribution. The ratio of the observed number of events to the expectation 
without oscillations is 0.74 ± 0.036 stclt ± 0.046 sys ±0.13t/j eor . 

The shape of the angular distribution of Fig. |]a has been tested with 
the hypothesis of no oscillation, giving a x 2 of 26.1 for 8 degrees of freedom. 
Assuming — > v T oscillations, the best \ 2 in the physical region of the 
oscillations parameters is 15.8 for Am 2 = 0.0025 eV 2 and sin 2 29 = 1. 

To test oscillation hypotheses, the independent probability for obtaining 
the number of events observed and the angular distribution for various pa- 
rameter values have been calculated. The value of Am 2 suggested from the 
shape of the angular distribution is similar to the value needed to obtain the 
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Figure 5: MACRO data. Measured and expected number of low energy events 
versus zenith angle; left: up stopping plus down semicontained; right: up semicon- 
tained. The solid lines are the predictions without oscillations; the dashed lines 
are the predictions assuming neutrino oscillations with the parameters suggested 
by the Up throughgoing sample. 



observed reduction in the number of events in the hypothesis of maximum 
mixing. Fig. ^b shows probability contours for oscillation parameters using 
the combination of probability for the number of events and the \ 2 °f the 
angular distribution. The maximum probability is 17%. The probability for 
no oscillations is 0.1%. 

The MC simulation for the low energy data uses the Bartol neutrino flux 
and the neutrino low energy cross sections of refBJ. The number of events 
and the angular distributions are compared with the predictions in Table ||] 
and Fig. [5[ The low energy data show a uniform deficit on the measured 
number of events over the whole angular distribution with respect to the 
predictions; there is good agreement with the predictions based on neutrino 
oscillations using the parameters obtained from the up throughgoing muon 
sample. 

6. Results from SuperKamiokande and Kamiokande 

SuperKamiokande is a large cylindrical water Cherenkov detector of 39 m 
diameter and 41 m height containing 50000 m 3 of water. The inner detector 
of 22500 m 3 is seen by 11146, 50-cm-diameter inner-facing phototubes. The 
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Figure 6: SuperKamiokande data. Zenith angle distributions of jU-like and e- 
like events for sub-GeV and multi-GeV data sets. Upward- going particles have 
cos0 < 0. The hatched regions are the Monte Carlo expectation for no oscillations 
normalized to the detector live-time. The solid lines are the best-fit expectations 
for — > v T oscillations with the overall flux normalization fitted as free parameter. 



2m thick outer layer of water acts as an anticoincidence and is seen by 1885 
smaller outward-facing photomultipliers. The ultra pure water has a light 
attenuation of almost 100 m. The detector is located in the Kamioka mine, 
Japan, under 2700 m.w.e.. For more details see the lecture by M. KoshibaEl. 

Because of the water index of reflection of 1.33, a relativistic charged 
particle (as the muon generated by a neutrino CC interaction) generates a 
42° forward cone of light. Instead electrons, being much lighter than muons, 
suffer electromagnetic showering and multiple scattering; therefore the ring 
of light is not so well denned as for the muons. The detector can distinguish 
single muons from single electrons with about 98% efficiency. Multiple ring 
events have been only partially used. The light intensity on the phototubes 
is a measure of the particle energy, and its direction is determined by the 
spatial and temporal pattern of phototube hits. 

The data presented at Neutrino'98 correspond to 33 kt ■ yr of data. This 
exposure collected 4353 fully contained (FC) events and 301 partially con- 
tained (PC) events. The FC events were separated into "sub-GeV" (with 
E V i S < 1.3 GeV) and "multi-GeV" (E vis > 1.33 GeV) samples; E vis is de- 
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Figure 7: SuperKamiokande data: (a) The ratio of the number of FC events 
to Monte Carlo events versus reconstructed L/E v . The points are the ratios of 
observed data to MC expectation in the absence of oscillations. The lower dashed 
line is the expected shape for Uu — > v T for Am 2 = 2.2 x 10~ 3 and sin 2 26 = 1. (b) 
The 68%, 90% and 99% C.L. contours for sin 2 26 and Am 2 for — > ^ T oscillations. 
The 90% C.L. contour obtained by the Kamiokande experiment is also shown. 



fined to be the energy of an electron that would produce the observed amount 
of Cherenkov light. E vis = 1.33 GeV corresponds to ~ 1.4 GeV/c. 

In a full-detector Monte Carlo simulation, 88% (96%) of sub-GeV e-like 
(//-like) events were u e (z/ M ) charged-current interactions and 84% (99%) of 
the multi-GeV e-like (//-like) events were v e [v^] CC interactions. PC events 
were estimated to be 98% charged-current interactions; hence, all the PC 
events were classified as //-like, and no single-ring requirement was made. 

The zenith angle distributions for e-like and //-like FC and PC events are 
shown in Fig. []. Note the near agreement of the e-like measured events with 
the MC predictions without oscillations, and instead the deviations from the 
no oscillation predictions for the //-like events. Significantly small values of 
the double ratio R = {^/e)data/{^/ e )MC in both the sub-GeV and multi-GeV 
samples were obtained. Several sources of systematic uncertainties in these 
measurement have been considered, but none of these can alter significantly 
the results. Moreover, the //-like data exhibit a strong up/down asymmetry 
in zenith angle (B) while no significant asymmetry is observed in the e-like 
data. The up/down ratio is expected to be near unity, really independent 
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of the flux model for E u > 1 GeV (above this value the effects due to the 
Earth' magnetic field on cosmic rays are small). The experimentally mea- 
sured up/down ratio is 0.54±g;g| ± 0.01 for the multi-GeV FC+PC /i-like 
events. 

SuperKamiokande estimated the oscillation parameters considering the R 
measurements and the zenith angle shapes separately. The 90% CL allowed 
regions for each case overlapped at 1 x 10 -3 < Am 2 < 4 x 10" 3 eV 2 for 
sin 2 26 = 1. Fig. [7|a shows the ratio of FC data to Monte Carlo for e-like 
and /i-like events with p > 400 MeV/c as a function of L/E u , compared to 
the expectation for — > v T oscillations with the best fit parameters from 
R versus zenith angle. While the e-like data show no significant variation 
in L/E u , the /i-like events show a significant deficit at large L/E v . At large 
L/E u , the have presumably undergone numerous oscillations and have 
averaged out to roughly half the initial rate. 

The final values for — » v T oscillations are sin 2 28 flT = 1 and Am 2 = 
0.0022 eV 2 . The contour plots for the neutrino oscillation parameters for 
Kamiokande and SuperKamiokande are shown in Fig. [?|b. Superkamiokande 
reported also data on up throughgoing muons, which agree with the predic- 
tions with an oscillated flux with the above parameters^. 

7. Conclusions 

A wealth of new data on atmospheric neutrinos was presented at the Neu- 
trino'98 Conference. The zenith angle distributions of atmospheric neutrino 
induced muons differ in shape and in absolute value from the ones predicted 
in the absence on neutrino oscillations. In the vertical upgoing direction 
there are about 50% deficits for low and high energy muon events. For v e 
induced electrons there is no strong deviation from prediction. The ratio of 
muons to electrons normalized to the respective MC predictions enhances 
the anomaly. 

The new data are in agreement with the hypothesis of two flavour v T 
oscillations, with maximum mixing and Am 2 ~ 0.0023 eV 2 ; the uncertainty 
in the Am 2 is relatively large. The 90% CL contours of MACRO and Su- 
perKamiokande overlap closely. 

The present experiments on atmospheric neutrinos are disappearance ex- 
periments; it would be nice to have a cross check with an appearance exper- 
iment . 
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One cannot exclude — > v s oscillations into a sterile v 3 . Answers to 
some detailed questions, like the shape of the high energy muon angular 
distribution may need more data. 

We would like to acknowledge the cooperation of many members of the 
MACRO collaboration, in particular of all the members of the neutrino work- 
ing group. 
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